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Abstract. In this paper we study the role of dynami- 
cal friction on the evolution of a population of large ob- 
jects (m > 10 22 g) at heliocentric distances > 70 AU 
in the Kuiper Belt. We show that the already flat dis- 
tribution of these objects must flatten further due to non- 
spherically symmetric distribution of matter in the Kuiper 
Belt. Moreover the dynamical drag, produced by dynam- 
ical friction, causes objects of masses > 10 24 g to lose an- 
gular momentum and to fall through more central regions 
in a timescale « 10 9 yr. This mechanism is able to trans- 
port inwards objects of the size of Pluto, supposing it was 
created beyond 50AU, according to a Stern & Colwell's 
(1997b) suggestion. 
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1. Introduction 

The current model for comets in the solar system sup- 
poses that a vast cloud of cometary objects orbits the 
Sun. This cloud consists of three components. The inner 
one, referred to as the Kuiper Belt (hereafter KB) (Edge- 
worth 1949; Kuiper 1951), is a disc like structure of > 10 10 
comets extending from 40 — 10 3 AU from the Sun (Weiss- 
man 1995; Luu et al. 1997). The KB has been proposed 
as the source of the Jupiter-family short-period (hereafter 
SP) comets. The second component, referred to as the 
Oort inner cloud, or the Hills cloud (Hills 1981), is sup- 
posed to be a disc, thicker than KB, containing 10 12 — 10 13 
objects lying - 10 3 - 2 x 10 4 AU fr om the Sun. It is pro- 
posed as a source of long-period (hereafter LP) and Halley- 
type SP comets (Levison 1996). The last component, the 
Oort cloud (Oort 1950), is a spherical cloud of 10 11 - 10 12 
cometary objects with nearly isotropic velocity distribu- 
tion extending from 2 x 10 4 to 2 x 10 5 AU. Even if the 
Oort cloud was considered in the past as the fundamental 
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reservoir of LP comets which have been brought into the 
inner solar system by perturbations due to the galactic 
tidal field, molecular clouds and passing stars, nowadays 
it has been shown that it can contribute only for a small 
part to the LP population of comets (Duncan et al. 1988; 
Wiegert & Tremaine 1997). 

Observational confirmation of the KB was first achieved 
with the discovery of object 1992QB1 by Jewitt & Luu 
(1993). To date over 40 KB objects (hereafter KBO) with 
diameters between 100 and 400 km have been discovered 
and the detection statistics obtained suggest that a com- 
plete ecliptic survey would reveal 7 x 10 4 such bodies 
orbiting between 30 and 50AU. Such a belt of distant 
icy planetesimals could be a more efficient source of SP 
comets than the Oort cloud (Fernandez 1980; Duncan 
et al. 1988). Dynamical simulations have shown that a 
cometary source with low initial inclination distribution 
was more consistent with the observed orbits of SP comets 
than the randomly distributed inclinations typical of the 
comets in the Oort cloud ( Quinn et al. 1990; Levison & 
Duncan 1993). According to other simulations the great- 
est part of objects in the KB should be stable for the age 
of the solar system. However if the population of the KB is 
~ 10 10 objects, weak gravitational perturbations provide 
a large enough influx to explain the current population 
of SP comets (Levison & Duncan 1993; Holman & Wis- 
dom 1993; Duncan et al. 1995). In particular Levison & 
Duncan (1993) and Holman & Wisdom (1993), studied 
the long term stability of test particles in low-eccentricity 
and low-inclination orbits beyond Neptune, subject only 
to the gravitational perturbations of the giant planets. 
They found orbital instability on timescales < 10 7 yr inte- 
rior to 33 — 34AU, regions of stability and instability in 
the range 34 — 43AU and stable orbits beyond 43AU. 
A study by Malhotra (1995a) showed that the KB is char- 
acterized by a highly non uniform distribution: most of 
the small bodies in the region between Neptune and 50AU 
would have been swept into narrow regions of orbital res- 
onance with Neptune (the 3:2 and 2:1 orbital resonances, 
respectively located at distances from the Sun of 39.4AU 
and 47.8AU). The orbital inclinations i of many of these 
objects would remain low (i < 10°) but the eccentric- 
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ities e would have values from 0.1 to 0.3. At the same 
time many of the trans-neptunian objects discovered lie 
in low-inclination orbits, as predicted by the dynamical 
models of Holman & Wisdom (1993) and Levison & Dun- 
can (1993). A more detailed analysis of this distribution 
reveals that most objects inside 42AU reside in highcr-e, i 
orbits locked in mean motion resonance with Neptune, but 
most objects beyond this distance reside in non-resonant 
orbits with significantly lower eccentricities and inclina- 
tions. 

After the previously quoted discovery of 100—200 km sized 
objects (Jewitt & Luu 1993; Jewitt & Luu 1995; Weiss- 
mann & Levison 1997), proving that the KB is populated, 
Cochran et al. (1995) have reported Hubble Space Tele- 
scope results giving the first direct evidence for comets in 
the KB. Cochran's observations imply that there is a large 
population (> 10 8 ) of Halley-sized objects (radii ~ 10 km) 
within <~ 40AU of the Sun, made up of low inclination ob- 
jects (i < 10°). At the time of the publication, Cochran's 
et al. (1995) results were criticized on two grounds: 

1) the detections were statistical in nature, and the au- 
thors were not able to fit orbits to their objects; 

2) the number of detections did not agree with extrapo- 
lation of the size distribution of large KBOs determined 
from early ground-based observations (but Weissman & 
Levison 1997 showed that Cochran's et al. 1995 results 
were in agreement with the number of KBOs needed to 
populate the Jupiter- family comets); moreover Brown et 
al. 1997 contended that detections reported in Cochran et 
al. (1995) were not possible, based on an analysis of the 
noise properties of the data. In a recent paper, Cochran 
et al. 1998, confirmed the early results by means of a new 
analysis. 

The spatial dimensions and mass distribution in the KB 
are poorly known. Yamamoto et al. (1994) have applied 
a planctcsimal model to the trans-neptunian region, find- 
ing that the maximum number density of the planetesimal 
population should be about at 100 — 200 AU and the plan- 
etesimal disc itself can extend up to distances « 10 3 AU. 
This is in agreement with detection by IRAS of discs 
around main sequence stars, Vega (Aumann et al. 1984), 
f3 Pictoris (Smith & Terrile 1984), extending to several 
hundred AU. From the available radio and infrared data, 
Beckwith & Sargent (1993) conclude that disc masses may 
range between 1O _3 M to 1M Q and extend from a few 
hundred AU to more than 10 3 AU. In short both theoret- 
ical arguments and observations strengthen the view that 
our solar system is surrounded by a flattened structure of 
planetesimals, extending perhaps to several hundred AU. 
Although originally it had been thought that the popula- 
tion might be collisionless, recent work (Stern 1995) has 
shown that the collisional effects cannot be neglected over 
4.5 Gyr. As shown by Stern (1995, 1996a,b) and Stern 
& Colwell (1997a, b), the collisional evolution is an im- 
portant evolutionary process in the disc as a whole, and 
moreover, it is likely to be the dominant evolutionary pro- 



cess beyond 42AU. In the case of larger planetesimals the 
evolution is connected to the energy loss due to dynami- 
cal friction, which transfers kinetic energy from the larger 
planetesimals to the smaller ones. This mechanism, in the 
early solar system, provides an energy source for the small 
planetesimals that is comparable to that provided by the 
viscous stirring process (Stewart & Wetherill 1988; Wei- 
denschilling et al. 1997). 

The objective of this paper is to examine the role of dy- 
namical friction, in the primordial KB, in the orbital evo- 
lution of the largest planetesimals that lie at a heliocentric 
distance > 70AU (at this distance the effects of the plan- 
ets decline rapidly to zero and only a small fraction of ob- 
jects is influenced by planetary perturbations - Wiegert & 
Trcmainc 1999; Stern & Colwell 1997b). While the impor- 
tance of dynamical friction in planetesimal dynamics was 
demonstrated in several papers, (Stewart & Kaula 1980; 
Horedt 1985; Stewart & Wetherill 1988) and in particular 
in the case of the planetary accumulation process, the role 
of this effect on the orbital evolution of the largest plan- 
etesimals and the consequent change of mass distribution 
in KB was never studied. 

The paper is organized as follows. In Sect. 2, we review 
the role of encounters and collisions in KB. In Sect. 3, we 
introduce the equations to calculate dynamical friction ef- 
fects. In Sect. 4 we describe how we use these equations to 
determine the evolution of the largest bodies population 
in KB. In Sect. 5 we discuss the results of the calculation 
and we also show (supposing the scenario proposed by 
Stern & Colwell 1997b of Pluto formation beyond 50 AU 
to be correct) how dynamical friction is able to transport 
an object of the size of Pluto from 50AU to the actual 
position. In Sect. 6 we give our conclusions. 

2. Encounters and collisions in KB 

As shown in several N-body simulations (Stern 1995), the 
structure of KB in the region 30 — 50AU is fundamentally 
due to two processes: 

— 1) dynamical erosion due to resonant interactions with 
Neptune (Holman & Wisdom 1993; Levison & Duncan 
1993; Duncan et al. 1995; Malhotra 1995a). In this re- 
gion the KB has a complex structure. Objects with 
perihelion distances ~ 35AU are unstable. For orbits 
with e > 0.1 and with semi-major axis a < 42 AU the 
only stable orbits are those in Neptune resonances. Be- 
tween 40AU and 42AU at low inclinations and between 
36AU and 39AU with % ~ 15° the orbits are unstable. 

— 2) collisional erosion (Stern 1995, 1996a,b; Stern & 
Colwell 1997a, b). Starting from a primordial disc hav- 
ing a mass of 4OM collisions are able to reduce its 
mass to 0.1 - O.3M in 10 9 yr if (e) > 0.1. 

Then the 30 — 50 AU zone is both collisionally and dynami- 
cally evolved, since dynamics acted to destabilize most or- 
bits with a < 42AU and were able to induce eccentricities 
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that caused collisions out to almost 50AU. The dynam- 
ically and collisionally evolved zone might extend as far 
as ~ 63AU, if Malhotra's (1995a) mean motion resonance 
sweeping mechanism is important. Beyond this region one 
expects there to be a collisionally evolved zone where ac- 
cretion has occurred but eccentricity perturbations by the 
giant planets have been too low to initiate erosion. Be- 
yond that region it is expected a primordial zone in which 
the accretion rates have hardly modified the initial pop- 
ulation of objects. Supposing that the last assumption is 
correct and that the radial distribution of heliocentric sur- 
face mass density in the disc, can be described by 
£(i?) oc R- 2 (Trcmaine 1990; Stern 1996a,b) and sup- 
posing that in the zone 30 — 50AU of the primordial disc 
40 — 50M e of matter was present, we expect ~ 35M® in 
the zone 70 - 100AU of the present KB. As the effects 
of planets are negligible for planetesimals beyond 70AU, 
the orbital motion there can be considered not far from 
Keplerian and circular (Brunini & Fernandez 1996). This 
last assumption is more strictly satisfied by the largest ob- 
jects, which should, through energy equipartition, evolve 
to the lowest eccentricity in the swarm (Stewart & Wether- 
ill 1988; Stern & Colwell 1997b) 

Although the main motion of KB objects (KBOs) is Ke- 
plerian rotation around the Sun, the motion is perturbed 
by encounters with other objects and by collisions. En- 
counters influence the structure of the system in several 
ways: 

a) Relaxation; 

b) Equipartition; 

c) Escape; 

d) Inelastic encounters. 

Each of the quoted effects has greater or smaller impor- 
tance in a system evolution according to the system char- 
acteristics. In the case of a system like KB, inelastic en- 
counters have a fundamental role because KBOs have on 
average a much smaller escape speed than the rms velocity 
dispersion, a. If r* is the radius of a KBO, = y / 2Gm/r* 

2 

is the escape speed from the KBO surface, O = is the 
Safronov number, n the number density and a the velocity 
dispersion, then the collision time t co \\ f° r a population of 
planetesimals with a Gaussian distribution in dispersion 
velocity is given by: 

(i) 



icon — 



lQ^narl{\ + 0) 



(Binncy & Tremaine 1987; Palmer et al. 1993). 
Within 1AU from the Sun, a population of a few hundred 
km-sized planetesimals, with several Earth masses in to- 
tal, would have i co ii — 10 4 yr, while at distances > 50AU 
t co \i becomes comparable with the age of the solar sys- 
tem. Indeed by means of N-body simulations, Stern (1995) 
showed how collisional evolution plays an important role 
through KB and that it has a dominant role at r > 42AU. 
The effect of collisions is that of inducing energy dissi- 
pation in the system but at the same time collisions are 



important in the growth of QB1 objects, Pluto-scale and 
larger objects starting from 1 to 10 km building blocks in 
a time that in some plausible circumstances is as little as 
~ 100 - 200Myr. 

Mutual gravitational scattering induces random velocity 
in two different ways: one is viscous stirring which con- 
verts solar gravitational energy into random kinetic en- 
ergy of planetesimals. Energy is transferred from circular, 
co-planar orbits with zero random velocities to eccentric, 
mutually inclined orbits with nonzero random velocities. 
The other is dynamical friction which transfers random ki- 
netic energy from the larger planetesimals to the smaller 
ones (Stewart & Wetherill 1988; Ida 1990; Ida & Makino 
1992; Palmer et al. 1993). Unlike viscous stirring, the ex- 
change of energy does not depend on the differential rota- 
tion of the mean flow for its existence. Dynamical friction 
would drive the system to a state of equipartition of kinetic 
energy but viscous stirring opposes this tendency. In the 
dispersive regime, the time scales of stirring and dynam- 
ical friction are almost equal to the two-body relaxation 
time 

1 
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r G is the gravitational radius, n the number density, 



G(ml + mi) 



(3) 



frmax being the largest impact parameter and |w| 2 the mean 
square velocity of the objects, which for the typical values 
of the parameters in the KB is of the order of w 10 11 and 
log A w 25. 

In conclusion the random velocity of the smaller planetesi- 
mals is increased by viscous stirring while the larger plan- 
etesimals suffer dynamical friction due to smaller plan- 
etesimals. 



3. Dynamical friction in KB 

The equation of motion of a KBO can be written as: 



'other (4) 



1 — F Q + Fpi ano t s + Ftido + F G CM + F sta rs + R + F, 

(Wiegert & Tremaine 1997). The term F Q represents the 
force per unit mass from the Sun, F p i anets that from plan- 
ets, F t ido that from the Galactic tide, Fqcm that from gi- 
ant molecular clouds, F sta rs that from passing stars, F ot hor 
that from other sources (e.g. non-gravitational forces), 
while R is the dissipative force (the sum of accretion and 
dynamical friction terms - sec Mclita & Woolfson 1996). If 
we consider KBOs at heliocentric distances > 70AU then 
Fpianots may be neglected. We also neglect the effects of 
non gravitational forces and the perturbations from Galac- 
tic tide, GCMs or stellar perturbations, which are impor- 
tant only for objects at heliocentric distances >> 100AU 
(Brunini & Fernandez 1996). We assume that the planetes- 
imals travel around the Sun in circular orbits and we study 
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the orbital evolution of KBOs after they reach a mass 
> 10 22 g. Moreover we suppose that the role of collisions 
for our KBOs at distances > 70AU can be neglected. We 
know that the role of collisions would be progressively less 
important with increasing distance from the Sun because 
the collision rate, nva' (a' is the collision cross section), 
decreases due to the decrease in the local space number 
density, n of KBOs and the local average crossing velocity, 
v, of the target body. As stated previously, using Eq. (|I]) 
at distances larger than 50 AU the collision time, f co ii, is of 
the order of the age of the solar system. Besides, the en- 
ergy damping is not dominated by collisional damping but 
by dynamical friction damping; also, artificially increasing 
the collisional damping the dynamics of the largest bodies 
is hardly changed (Kokubo & Ida 1998). 
To take account of dynamical friction we need a suitable 
formula for a disk- like structure such as KB. Following 
Chandrasekhar & von Neumann's (1942) method, the fric- 
tional force which is experienced by a body of mass mi, 
moving through a homogeneous and isotropic distribution 
of lighter particles of mass m,2 , having a velocity distribu- 
tion n(v2) is given by: 



-47777117712 (rTJl + m 2 )G 2 



Vi 

n(v 2 )dv 2 —T log A 



(5) 



(Chandrasekhar 1943); where log A is the Coulomb loga- 
rithm, 777! and 7772 are, respectively, the masses of the test 
particle and that of the field one, and v\ and t>2 the re- 
spective velocity, n(y2)dv2 is the number of field particles 
with velocities between v 2 , v 2 + dv 2 . 

If the velocity distribution is Maxwellian Eq. (I|) becomes: 



-A^nmximx + m 2 )pG 2 — i 
v i 



logA[er/(A) - 2X cxp(-A 2 )/ yf{v)} 



(6) 



(Chandrasekhar 1943, Binney & Tremaine 1987), where p 
is the density of field particles and X = v\j \] (2a), a 
being the velocity dispersion. Eq. (||) cannot be used for 
systems not spherically symmetric except for the case of 
objects moving in the equatorial plane of an axisymmetric 
distribution of matter. These objects, in fact, have no way 
of perceiving that the potential in which they move is not 
spherically symmetric. 

We know that KB is a disc and consequently for objects 
moving away from the disc plane we need a more gen- 
eral formula than Eq. (^). Moreover dynamical friction 
in discs differs from that in spherical isotropic three di- 
mensional systems. First, in a disc close encounters give 
a contribution to the friction that is comparable to that 
of distant encounters (Donner & Sundelius 1993; Palmer 
et al. 1993). Collective effects in a disc are much stronger 
than in a three-dimensional system. The velocity disper- 
sion of particles in a disc potential is anisotropic. N-body 
simulations and observations show that the radial compo- 
nent of the dispersion, ctr, and the vertical one, er z , are 



characterized by a ratio cr/<t z ~ 0.5 for planetesimals 
in a Keplerian disc (Ida et al. 1993). The velocity dis- 
persion evolves through gravitational scattering between 
particles. Gravitational scattering between particles trans- 
fers the energy of the systematic rotation to the random 
motion (Stewart & Wetherill 1988). In other words the ve- 
locity distribution of a Keplerian particle disc is ellipsoidal 
with ratio 2:1 between the radial and orthogonal (z) di- 
rections (Stewart & Wetherill 1988). According with what 
previously told, we assume that the matter-distribution is 
disc-shaped, having a velocity distribution: 



n(v, x) = n(x) 



3/2 



exp 



2a 



2a\ 



a 2 a± 



(7) 



(Hornung & al. 1985, Stewart & Wetherill 1988) where v\\ 
and crii are the velocity and the velocity dispersion in the 
direction parallel to the plane while v± and a± are the 
same in the perpendicular direction. We suppose that cry 
and a j_ are constants and their ratio is simply taken to be 
2:1. Then according to Chandrasekhar (1968) and Binney 
(1977) we may write the force components as: 



*ll = 



fc|l«i|l = B||Wi|| • 
2\f2~nnG 2 lg Ami 7772 (777 1 



777 2 



)v / r^ 2 



aja± 



F± = kxvi± = B±vi x ■ 

2\f2~nnG 2 \ghmim 2 (mi + 7772) \/\ — e 2 — j- 



a 2 a± 



(8) 



where 



B n = 



dqeup 



2a\ 1 + q 2a 2 ,l-e 2 + q 



(l + q) 2 (l-ei + q) 1/2 



(10) 



B 1 



dq exp 



i± 



2a 2 1 



2a\ 1 



e z + q 



1 



(l + q) 2 (l-e 2 + q) 



3/2 



and 



(11) 



(12) 



while 77 is the average spatial density. When B± > Bu the 
drag caused by dynamical friction will tend to increase the 
anisotropy of the velocity distribution of the test particles. 
The frictional drag on the test particles may be written: 




F = || ey - k±v 1 ±e A 



(13) 
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where e» and ej_ are two unity vectors parallel and per- 
pendicular to the disc plane. 

This result differs from the classical Chandrasekhar (1943) 
formula. Chandrasekhar's result tells that dynamical fric- 
tion force, F, is always directed as —v. This means that if 
a massive body moves, for example, in a disc in a plane dif- 
ferent from the symmetry plane, dynamical friction causes 
it to spiral through the center of the mass distribution al- 
ways remaning on its own plane. It shall reach the disc 
plane only when it reaches the centre of the distribution. 
Eq. ( [l3|) shows that a massive object suffers drags — k\\Vi« 
and — k±vi±_ in the directions within and perpendicular 
to the equatorial plane of the distribution. This means 
that the object shall find itself confined to the plane of 
the disc before it reaches the centre of the distribution 
(this means that also inclinations are damped). In other 
words the dynamical drag experienced by an object of 
mass mi moving through a non-spherical distribution of 
less massive objects of mass mi is not directed in the di- 
rection of the relative movement of the massive particle 
and the centre of mass of the less massive objects (as in 
the case of spherically symmetric distribution of matter). 
As a consequence the already flat distribution of more 
massive objects will be further flattened during the evo- 
lution of the system (Binney 1977). The objects lying in 
the plane, as previously told, have no way to perceive that 
they are moving in a not spherically symmetric potential. 
Hence we expect that the dynamical drag is directed in 
the direction opposite to the motion of the particle: 



(Binney & Tremaine 1987; Wiegert & Tremaine 1997): 



-fcii ('lien 



4. Parameters used in the simulation 



(14) 



To calculate the effects of dynamical friction, introduced 
in the previous section, on a disc-like structure as KB, we 
cannot use the classical Chandrasekhar's (1943) theory 
but we need equations specified for distributions like KB. 
These equations were written in the previous section (Eqs. 
H — [n]) . To calculate the effect of dynamical friction on the 
orbital evolution of the largest bodies we suppose that a± 
and an are constants and that a\\—2a±. We need also the 
mass density distribution in the disc. We assume a helio- 
centric (R) distribution in surface mass density S oc R~ 2 
and a total primordial mass, M = 50M^ in the region 
30 - 50AU (Stern & Colwell 1997b). To evaluate the dy- 
namical friction force we need the spatial distribution of 
the field objects, n. To reproduce the quoted surface den- 
sity we need a number density decreasing as a function of 
distance from the Sun, R, to the third power (Levison & 
Duncan 1990): 



n R 



(15) 



n(R, z) 



-in 



aR 2 + (a + 3Vz 2 + b 2 ) (a + ^z 2 + b 2 f 
R 2 + (a+ Vz 2 + b 2 ) 2 ] (z 2 + b 2 ) 3/2 



(16) 



We remember also that the KB is a disc and then we 
use the mass distribution given by a Myiamoto-Nagai disc 



which in the disc plane reproduces the quoted surface den- 
sity and Eq. @ [n(R,0) cx i?~ 3 ]. Here M is the disc 
mass, a and b are parameters describing the disc charac- 
teristic radius and thickness. Because there is presently 
no information on the way in which ensemble-averaged 
inclinations ((i)) and eccentricities ((e)) vary in the KB, 
we adopt a disc with a width (i) = |(e) (Stern 1996b). 
The equations of motion were integrated in heliocentric 
coordinates using the Bulirsch-Stoer method. 

5. Results 

The model described was integrated for several values of 
masses, starting from m = 10 22 g, supposing the KBOs 
move on circular orbits. We studied the motion of KBOs 
both on inclined orbits, in order to study the evolution 
of the inclination with time, and on orbits on the plane 
of the disc, in order to study the drift of the KBOs from 
their initial position. The masses of the KBOs, M, were 
considered constant during the whole integration in order 
to reduce the number of differential equations to solve. 
Moreover we use M << Nm and m « M, N and m 
being the total number and the mass of the swarm of field 
particles in which the KBO moves. The assumption that 
field particles have all equal masses, m, does not affect the 
results, since dynamical friction does not depend on the 
individual masses of these particles but on their overall 
density. The results of our calculations are shown in Figs. 
1 - 4. 

In Fig. 1 we plot the values of (z) / z a versus time for plan- 
etesimals having masses 10 22 g and 10 23 g. The z coordinate 
is orthogonal to the plane of the disc and z Q — 7AU while 
t a = 10 8 yr. The brackets () are mean values obtained av- 
eraging over a suitable number of orbital oscillations. As 
shown in the figure the decay of the inclination of the plan- 
etesimal of 10 22 g (dotted line) has a timescale larger than 
the age of the solar system and on the order of 3 x 10 10 yr. 
The inclination of the more massive planetesimal (full line) 
decays in a time ~ 3 x 10 9 yr. This is due to the fact that 
dynamical friction effects increase with the mass M of the 
KBO. In fact (see Fig. 2) when the mass of the KBO is 
10 25 g the decay time reduces to 3 x 10 8 yr. Dynamical fric- 
tion makes the orbits of KBOs undergo some collimation 
along the z direction, characterized by a low value of the 
dispersion velocity. This was expected because it is known 
that the larger the velocity dispersion along the direction 
of motion, the lesser the effect of dynamical friction (Pesce 
et al. 1992). Binney (1977) has found an efficient collima- 
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Fig. 1. Height from the plane versus time for a planetesi- 
mal of 10 22 g, dotted line, and one of 10 23 g, full line. Time 
is measured in units of t = 10 8 yr while height is measured 
in z Q = 7AU. 

tion of orbits along the main axis of the velocity dispersion 
tensor in the case of an anisotropic axisymmetric system, 
in which the principal velocity dispersions have constant 
values. 

The result obtained is in agreement with previous studies 
of the damping of the inclinations of very massive objects 
by Ida (1990), Ida & Makino (1992) and Melita & Woolf- 
son (1996). In fact in the semi-analytical theory by Ida 
(1990) the timescale for inclinations damping due to dy- 
namical friction is almost equal to the two-body relaxation 
time 



Fig. 2. Height from the plane versus time for a planetes- 
imal of 10 25 g. Time is measured in units of t = 10 8 yr 
while height is measured in z — 7AU. 

plane radial heliocentric distance of the planctesimal while 
r Q = 70AU and t Q — 10 10 yr. As shown, the time required 
to a planetesimal of the quoted mass to reach 40AU is ~ 
3 x 10 10 yr, which is an order of magnitude larger than the 
damping timescale. This is due, in agreement with what 
previously told, to the fact that in the plane the dispersion 
velocity is larger than that in the z direction. Increasing 
the mass to 10 25 g the time needed for a planetesimal to 
reach 40AU decreases to 4 x 10 8 yr (here t = 10 6 yr) (see 
Fig. 4). The threshold planetesimal mass, Mthreshoidi that 
starts orbital migration is ~ 10 24 g. This mass scales with 
the disk density as: 



1 



1 



damp - T 2B ~ ^ 2M2 lnA/ ^ ~ SGW fi , nA M threshold 



10 ; 



24 



(p/TO-^g/cm 3 ) 



(18) 



where S is the surface density of small objects and fl is 
the Kcplcrian frequency The timescales given by Ida's 
theory for our planetesimals of 10 22 g, 10 23 g, 10 25 g are 
respectively 4 x 10 10 yr, 4 x 10 9 yr, 4 x 10 8 yr, in agreement 
with our results. 

The effect of dynamical friction on the semi-major axis is 
plotted in Fig. 3 and Fig. 4. In Fig. 3 we plot (r) jr Q versus 
time for a planetesimal of mass 10 23 g. Here r is the in- 



We recall that we do not take into account the effects of 
the planets because we considered planetesimals initially 
at distances > 70AU, but when the planetesimal moves 
towards the region of influence of planets the role of these 
must be taken into account. 

We have some difficulty to compare this result with previ- 
ous studies because the problem of the decay of the semi- 
major axis has not been particularly studied. So far, many 
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Fig. 3. Heliocentric distances in the disk plane for a plan- 
etesimal of 10 23 g. Time is measured in units of t — 10 10 yr 
while distances are measured in r = 70AU. 

people have assumed a priori that radial migration due to 
dynamical friction is much slower than damping of veloc- 
ity dispersion due to dynamical friction. Therefore most 
of studies of dynamical friction were concerned only with 
damping of velocity dispersion (damping of the eccentric- 
ity, e, and inclination, i), adopting local coordinates. An- 
alytical work by Stewart & Wetherill (1988) and by Ida 
(1990), adopted local coordinates. N-body simulation by 
Ida & Makino (1992) adopted non-local coordinates, but 
did not investigate radial migration. Only the density wave 
approach by Goldreich & Tremaine (1979, 1980) and Ward 
(1986) considered radial migration. However, the relation 
of this approach to particle orbit approach is not clear. 
Furthermore, a few numerical simulation has been devoted 
to investigate radial migration. 

In any case we shall compare our result with that by Ward 
(1986) supposing that this density approach describes cor- 
rectly the radial migration. Following Ward (1986) the 
characteristic orbital decay time of a disc perturbcr is 
given by: 

Tdcc= n\c\^z^ (^n) (19) 



Fig. 4. Heliocentric distances in the disk plane for a plan- 
etesimal of 10 25 g. Time is measured in units of t Q — 10 6 yr 
while distances are measured in r = 70AU. 

where fx = M/M@, and the nondimensional factor C, de- 
pending on the disc's surface density gradient, k, and the 
adiabatic index, s, for a disc with Q = oo is ~ —18, c is 
the gas sound speed (in a planetcsimals system this must 
be replaced by the velocity dispersion). The timescale for 
the perturber to drift out of a region of radius r is given 
by (Ward 1986): 

Tdrift — -^Trfec (20) 

Using Q — y/GM & /r 3 , calculating S supposing that the 
disc mass is uniformly spread in the region 40 — 70 AU and 
c = 20000 cm/s we find Tdrift — 7 x 10 8 yr for a planetesi- 
mal having M = 10 25 g, a little larger than the value pre- 
viously found. We have to remember that Ward's (1986) 
model supposes that the solar nebula is two-dimensional 
and that in a finite thickness disc other damping mech- 
anisms may come into play perhaps invalidating Ward's 
result. 

This last calculation shows that in a timescale less than 
the age of the solar system, objects of the mass of Pluto 
~ O.OO2M0 may move from the region > 50AU towards 
the position actually occupied by the planet. This opens 
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a third possibility to the standard scenarios for Pluto for- 
mation (in situ formation, formation at 20 — 30AU and 
transport outwards) namely that Pluto was created be- 
yond 50AU, and then transported inwards. One of the 
problems of Pluto formation in the first two scenarios is 
that the growth of Neptune caused accretion to be inhib- 
ited and then it is necessary that the time for accreting 
an object of Pluto mass was shorter than the timescale 
of Neptune formation. This problem is not present in the 
third scenario because at heliocentric distances larger than 
50AU Neptune never induced significant eccentricities on 
most orbits in the region. Hence the dynamical conditions 
necessary for growth may have persisted for the whole age 
of the solar system and consequently Pluto could have 
formed later than in the other two quoted scenarios. The 
mechanism responsible for the transport of Pluto from the 
quoted region to that nowadays occupied might have been 
dynamical friction. 

Three possible objections to this last model are: 

1) Moving to greater and greater distances both the disc 
density and the velocity dispersions decrease and conse- 
quently the accretion times increase. Can a Pluto-scale 
body form at distances of 70 AU? 

2) A possible explanation for Pluto's orbital parameters 
is connected to outwards migration of Neptune (Malho- 
tra 1993). If this works for Neptune, could it also work 
for Pluto-like objects in the KB? How much might Pluto 
have moved out? Could that compensate for the effect of 
the frictional motion? 

3) How can one explain the odd orbital parameters of 
Pluto's orbit (e = 0.25, i — 17 degrees to the ecliptic)? 
Surely, as made clear by the first objection, formation of 
large bodies is more and more difficult moving away from 
the inner parts of Solar system. Although growth times at 
70 AU are about 4-5 times longer than at 40 AU, Pluto- 
mass bodies can indeed be grown at this distance, from 
1 to 10 km building blocks, in ~ lGyr, if the mean disc 
eccentricity, < e >~ 0.001, and if the KB mass interior to 
50 AU was, as previously stated in Sect. 4, 30-50 and 
continued outwards with S cx R~ 2 (see Stern & Colwell 
1997b). 

The answer to the second objection is the following. 
Orbital migration of a planet can be accomplished through 
two mechanisms. In the first mechanism a planet and the 
circumstellar disc interact tidally which results in angular 
momentum transfer between the disc and the planet (e.g. 
Goldreich & Trcmainc 1980; Ward 1997). The planet's mo- 
tion in the disc excites density waves both interior and ex- 
terior to the planet. If the planet is large enough (at least 
several Earth masses), it is able to open and sustain a gap. 
It establishes a barrier to any radial disc flow due to vis- 
cous diffusion and it becomes locked to the disc and must 
ultimately share its fate (this is known as type II drift). 
In this case both inwards and outwards planet migration 
are allowed. In fact in a viscous disc, gas inside the ra- 
dius of maximum viscous stress, r mvs , drifts inwards as it 



loses angular momentum while gas outside r mvs expands 
outwards as it receives angular momentum (Lyndcn-Bell 
& Pringle 1974). Neptune's outwards migration is due to 
the fact that the gas in the Neptune forming region has a 
tendency to migrate outwards (Ruden & Lin 1986). 
If the planet is not able to sustain a gap, the net torque 
from the disc is still not zero and it migrates inwards in a 
shorter timescale (type I drift). 

Pluto being a low mass planet it can migrate only by 
means of type I drift, this means that it can only migrate 
inwards. 

In the second mechanism a planet can undergo orbital 
migration as a consequence of gravitational scattering be- 
tween itself and residuals planetesimals. If a planctcsimal 
in a near-circular orbit similar to that of the planet is 
ejected into a Solar system escape orbit, the planet suffers 
a loss of orbital angular momentum and a corresponding 
change of orbital radius. Conversely, planetesimals scat- 
tered inwards would cause an increase of orbital radius 
and angular momentum of the planet. A single massive 
planet scattering a population of planetesimals in near- 
circular orbits in the vicinity of its own orbit would suf- 
fer no net change of orbital radius as it scatters approxi- 
mately equal numbers of planetesimals inwards and out- 
wards. However in some peculiar situations, such as that 
encountered in the region of Jovian planets, things go dif- 
ferently from this picture (Fernandez & Ip 1984). In par- 
ticular, as Jupiter preferentially removes the inwards scat- 
tered Neptune planetesimals, the planctcsimal population 
encountering Neptune at later times is increasingly biased 
towards objects with specific angular momentum larger 
than Neptune's. Encounters with this planctcsimal pop- 
ulation produce a net gain of angular momentum, hence 
an increase in its orbital radius. Evidently this situation 
is not the one present in the outer Solar system region, 
occupied by Pluto in our model. In other words, there is 
no reason to suppose that Pluto moved outwards like Nep- 
tune. 

For what concerns the third objection, a possible answer 
is that Pluto gained high eccentricity and inclination in a 
similar way to that described by Malhotra (1993, 1995b). 
There has, of course, been much speculation as to the ori- 
gin of the extraordinary orbit of Pluto (Lyttleton 1936; 
Farinella et al. 1979; Olsson-Steel 1988; Malhotra 1993). 
All but one (Malhotra (1993, 1995b)) of these specula- 
tions require one or more low-probability "catastrophic" 
events. In Malhotra's (1993, 1995b) model, Neptune's or- 
bit may have expanded considerably, and its exterior or- 
bital resonances would have swept through a large region 
of trans-Neptunian space. During these resonances sweep- 
ing, Pluto could have been captured into the 3:2 orbital 
period resonance with Neptune and its eccentricity and 
inclination would have been pumped up during the sub- 
sequent evolution. 

The phenomenon of capture into resonance as result of 
some dissipative forces is common in nature. Weiden- 
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schilling & Davies (1985) studied resonance trapping of 
planetesimals by a protoplanet in association with gas 
drag. Many of the characteristics of this effect have been 
studied. Patterson (1987) and Beauge et al. (1994) have 
investigated its cosmogonic implications. The stability of 
the orbits and capture probabilities have been studied by 
Beauge & Ferraz-Mello (1993), Gomes (1995). Melita & 
Woolfson (1996) showed that a three-body system (Sun 
and two planets) under the influence of both accretion and 
dynamical friction forces, evolve into planetary resonance 
when the inner body is more massive. In general capture 
into a stable orbit-orbit resonance is possible when the or- 
bits of two bodies approach each other as a result of the 
action of some dissipative process. The transition from a 
non-resonant to a resonant orbit depends sensitively upon 
initial conditions and the rate of orbital evolution due 
to the dissipative effects. Borderies & Goldreich (1984) 
showed that for single resonance and in the limit of slow 
" adiabatic" orbit evolution, the probability of capture for 
the 3:2 Neptune resonance is 100% for initial eccentricity 
less than ~ 0.03 and reduces to less than 10% for initial 
eccentricities exceeding 0.15 (see Malhotra 1995b). 
In our model Pluto reaches its actual position in a time (> 
10 8 yr) larger than Neptune's orbital migration timescale 
(10 6 - 10 7 yr; see Ida et al. 1999). Then Neptune was 
in its actual position when Pluto reached its own. Pluto, 
before the resonance encounter, Pluto has an initial low 
eccentricity ~ 0.03, because of eccentricity and inclina- 
tion damping due to dynamical friction, and its migration 
is very slow (> 10 8 yr). The capture probability is then 
very high. The increase in eccentricity and inclination is 
naturally explained by Malhotra's theory. 



6. Conclusions 

In this paper we studied how dynamical friction due to 
small planetesimals influence the evolution of KBOs hav- 
ing masses larger than 10 22 g. We find that mean eccentric- 
ity of large mass particles is reduced by dynamical friction 
by small mass particles in timescales shorter than the age 
of the solar system for objects of mass equal or larger than 
10 23 g. We also studied the effect of dynamical friction on 
the evolution of the semi-major axis of the largest plan- 
etesimals. We find that even if dynamical friction is less 
effective in transferring planetesimals towards the inner 
part of the solar system, with respect to the damping of 
inclinations that it is able to produce, the timescale for 
radial migration is shorter than the age of the solar sys- 
tem for large enough masses (> 10 24 g). 
Finally, our calculation show the dynamical friction may 
be the mechanism responsible for the transport of objects 
like Pluto from regions with r > 50AU towards the posi- 
tion nowdays occupied and this opens a third possibility 
for Pluto formation that eliminates the problem of the 
Neptune formation. 
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